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Abstract 

We derive the first terms in the amphtude of lepton pair production in the 
Coulomb fields of two relativistic heavy ions. Using the Sudakov technique, 
which very simplify the calculations in momentum space for the processes at high 
energies, we get the compact analytical expressions for differential cross section 
of the process under consideration in the lowest order in fine structure constant 
(Born approximation) valid for any momentum transfer and in a wide kinematics 
region for produced particles. Exploiting the same technique we consider the 
next terms of perturbation series (up to fourth order in fine structure constant) 
and investigate their energy dependence and limiting cases. It has been shown 
that taking in account all relevant terms in corresponding order one obtains the 
expressions which are gauge invariant and finite. We estimate the contribution of 
the Coulomb corrections to the total cross section and discuss the cancellations 
of the different terms which holds in the total cross section. 



1 Introduction 

The process of lepton pair production in collision of relativistic heavy ions attracted 
a rising interest in the past years which is connected mainly with operation of RHIC 
( Lorentz factor j = E/M = 100) and LHC (7 = 3000). The total cross section of the 
production of lepton pair in the collision of two relativistic particles 

A + B^e+e'+A + B (1) 
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in its lowest order in fine structure constant = |^ = 3^ has been known long ago 
^ 0. Over the last decade a lot of works have been done to investigate this process for 
the case when A and B are heavy ions (see e.g. |Q and references there in). When heavy 
ions collide at relativistic velocities their electromagnetic fields are Lorentz contracted 
and the process (|l]) can be viewed as pair production in two external fields which can 
be represented by classical electromagnetic potentials. Such approach allows 0, || to 
investigate the impact parameter dependence of this process and calculate its cross 
section with higher accuracy, than predicted by widely exploited Weizsacker- Williams 
(WW) approximation. This works are based on the Born approximation (see Fig. [l|) 
but even in this case the obtained expressions are complex, so that the certain "state 
of art" efforts in their computing are required 0. 

The topical problem is to take into account the Coulomb corrections (CC) to the 
amplitude of this process (see Fig. ^, |^). In spite of the fact that CC are determined 
by the powers of fine structure constant a, it enters the amplitude in combination 
(aZi)"^^ {aZ2)"''^ {aZiZ2)^^ , where Zi, Z2 are the colliding ions charge numbers, ni+n2 > 
2 is the number of exchanged photons between nuclei and the produced pair, ns > is 
the number of photons exchanged between the colliding nuclei. For heavy ion collisions 
the parameter aZ is of the order of unity and the problem arise to take into account the 
effect of CC. The exact knowledge of total yield of lepton pairs in heavy ion collisions is 
important issue, because the pair production consist a huge background in experiments 
with relativistic heavy ions. For example the total cross section for the process (1) at 
RHIC energies is tens kilobarns for the heavy ion collisions. At LHC energies this 
quantity becomes a hundreds kilobarns according to the Racah formula [0 . Moreover, 
the intensive pair production can destroy the ion beams circulating in the accelerator 
as a result of electron capture by heavy ion (see e.g. 0). 

Recently the problem of CC in the process (1) was investigated in the number of 
papers [§, ^ |10|, |TT[ with the unexpected result. The authors |10| solved the Dirac 
equation for an electron in the Coulomb field of two highly relativistic nuclei. Having 
used the crossing symmetry property, they connected the amplitude for electron scat- 
tering in the Coulomb fields of two heavy ions with amplitude for pair production^. 
As a result the total cross section of the process (1) turn out to coincide with its Born 
approximation, so that CC haven't impact on the total lepton pairs yield. The same 
result have been obtained in the work |TT]] in which the authors use the eikonal approx- 
imation for electron scattering at high energies. They took advantage of the common 
wisdom that the interaction of lepton pair with Coulomb fields of two highly relativistic 
nucleus can be represented as a product of the eikonal amplitudes for interaction of 
electron and positron with a nucleus A and B separately. Such approach allows one 
to take into account the contribution of CC in amphtude of the process (^ with the 
quoted above result: the total cross section exactly coincide with the Born term. This 
works induced the series of critical papers O, where it was explained why this 



result is so unexpected and how assumptions which have been done through its deriva- 
tion can affect this result. As was observed in [O the main contribution to the CC 
in the process (|l|) comes from the case when one of the colliding ions radiate a single 
photon and the produced pair interact with another nucleus by means of the arbitrary 
number of photons. In fact the contributions of such type can be obtained with the 



^The similar trick was exploited many years ago [l^ to get the probability of pair production in 
the condensed media from the probability of the electron bremsstrahlmig in the matter. 
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use of common WW approximation, which leads to appearance of large logarithms in 
energy for such terms, which can't be cancelled by the next order of CC where such 
logarithms are absent. On the other side it has been shown that the change in the 
order of integration in regularized integrals (the Coulomb phase for electron scattering 



diverges, so it needs regularization) completely change the final result. In |]10| the two 
approaches employed in investigation of this problem were analyzed in detail, with the 
conclusion that the perturbation approach in fact correspond to the exclusive process of 
the coherent production of lepton pairs in the collisions of relativistic ions (i.e. reaction 
(|l])). As to the approach based on retarded solutions of the Dirac equation it gives the 
inclusive cross section for production of arbitrary number of the lepton pairs, i.e. in 
this way one obtain the contribution of inclusive pairs production. The reason of such 
unexpected result from our point of view is the inapplicability of crossing symmetry 
property to the higher order (beyond the Born approximation) terms, the fact, which 
is known long ago (see e.g. [|1^) Thus the issue of calculation of the CC in the correct 
way in the process (|T]) is an open task up to now and it turn out to be more complex 
problem that it seems from the first glance. All this stimulated us to investigate the 
CC contribution in the perturbative expansion of amplitude for exclusive process under 
discussion in more details, with the aim to understand how one can take into account 
the CC in heavy ion collisions in the framework of perturbative QED. 

Our paper is organized as follows. In the next section 2 we present the amplitude 
for the lepton pair creation in the collision of relativistic particles in the Born approx- 
imation. We are restricted by kinematics of pionization region which gives the main 
contribution to the cross section, so that the pair production doesn't influence the 
kinematics of the colliding particles. Moreover, we supposed that momenta transfer 
between the produced pair and initial particles are much smaller in comparison with 
the energy of initial particles, e.g. we consider the case of their quasi-elastic scattering. 
We use the powerful Sudakov technique for the evaluation of the Feynman diagrams in 
momentum representation which allows to obtain the results with power accuracy in 
energy. The compact analytical expressions have been obtained for the Born amplitude 
and differential cross section of the process (1) which allows one to calculate it in the 
wide interval of kinematic variables. In the limiting case of quasi-real photons the well 
known results of WW approximation can be easily obtained from this expressions. 

In the section 3 we consider all possibilities for pair production by three photon 
exchange, so called amplitudes of the type M^^y We present the general formula for its 
contribution to the singlet triplet pair production amplitude and relevant cross section. 
Our result appears to be the generalization of the known result for the amplitude of 
the process of pair photoproduction in the odd charge state derived in DTG). We have 
done the numerical estimate for the relevant contribution to the total cross section in 
the main (double logarithmic) approximation. 

In the section 4 we present the results of the amplitude computation and corre- 
sponding contributions in the cross section from processes of two photon exchange 
between each ion and the created pair (fourth order terms in fine structure constant) 
and the case, when one of the ions is attached to the pair by single photon and another 
ion is connected with a pair by three photons in all possible ways. In the last case, our 
result appears to be generalization of the known result in the case of photoproduction 



of the pair in even charge state [0. We show that in any order the obtained expres- 
sions are gauge invariant and finite. We analyze the main contributions from obtained 
CC to the total cross section and show that the number of remarkable cancellation 
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between the different CC and their interference with the Born term take place. 

In the section 5 we shortly listed the main results obtained in the paper. 

In the Appendix A we give the results of analytical integration on loop momenta 
for spin-singlet and spin-triplet parts of CC, which allows one to obtain the first CC 
in the analytical form. 

In the Appendix B we present the spin structures which determines the fourth order 
terms in the amplitude of the process under consideration. 



2 The Born approximation 

We are interested in the exclusive coherent process of lepton pair production in the 
collision of two relativistic nuclei with charge numbers Zx , Z2 

A{vx) + B{j>2) ^ e+(g+) + e-(g_) + + B{j^^). (2) 

We define the usual c.m.s. total energy of colliding nuclei s = {pi +^2)^ = 471727711/712 
where rrii and 7i = ^ are the masses and Lorentz factors of the colliding nuclei. 
Keeping in mind the fast decrease of 7*7* —>■ e~^e~ cross section with the momentum 
transfer, i.e. the photons virtuality and with the invariant mass of the pair Sp = 
(g+ + we will work in the kinematics 

s = (p^+p^y > \ql\^ |g2|^ Sp. 



o2 



Later on we calculate the amplitudes neglecting the pieces of the order so our re- 
sults are valid for relativistic particles with power accuracy. The total cross section of 
the reaction in the lowest order in fine structure constant (Born approximation) is 
known long time ago 0. Nevertheless we carry out the derivation of amplitude and 
differential cross section in this approximation using the Sudakov technique, which is 
very useful for calculation of Feynman diagrams (FD) for the processes at high ener- 
gies. This allow us to illustrate the method and approximations which we will use later 
on in the CC derivation. Moreover we will obtain the compact analytical expressions, 
which are valid not only for the heavy ion collisions, but for the lepton pair production 
in the Coulomb fields of any relativistic particles (hadrons, leptons). 

The lowest order FD describing the process (|2|) are depicted in the Fig. |1]. The 
corresponding amplitude read 

(1) ^ _. (47ra) ^^^^ u{p[)-f^u{pi)u{p2)-fMP2)gf,aguf3u{9-)Taf3v{q+) (3) 
qi-q+ + m q2 - q+ + m 

J-af} = 7/37 ^ + 7a 7 ^ ^7/3- 

Using the Gribov's decomposition of metric tensor g^j^^, = g^^ + ^{PiP2 +P1P2) it is easy 
to show that at high energies the main contribution is given by its longitudinal part 
g^iaQufi (7) PiP2PiP2- As a result the expression can be cast in the form 

M« = -^s^^^^4^H^-)RU-(^.)N^N2, (4) 



11 1 

^(J) = -PlPiT^v: Ni = -U{p[)p2u{pi), N2 = -U{P2)P1U{P2). 
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For spinless projectile and target as well as for the colliding fermions at definite chirality 
state \Ni\ = \N2\ = 1, so that later onQ we will set iVi = = 1. It should be noted 
that if we consider the heavy ion's dimensions it is enough to multiply the amplitude 
(^) by the colliding ion's formfactors F{qf). 

In what follows we will use the standard Sudakov expansion of all the momenta in 
the two almost light-light vectors pi, p2 and the two-dimensional transverse component 
q±Pi = q±P2 = 

qi = aip2 + Pipi + qi±, q2 = 012P2 + P2P1 + q2±, q± = (y±p2 + P±pi + q±±, (5) 

Pi = Pi - P2 — , P2=P2-P1 — , Pl=P2 = O {^] , 2piP2 ^ 2piP2- 

s s \ s J 

From the on-mass shell condition, for instance p'l^ = {pi — qiY = rnf, p'2^ = (p2 — 92)^ = 
ml, one can easily see that /3i ~ 0:2 ~ ai ~ ~ ^ which allows to neglect, where 
it is possible, the values of ai, (32 which greatly simplify calculations and allows to 
obtain the results with the power accuracy. In this approximation we have 

gi~/5ipi + gi±, g2~«2P2 + g2±, (6) 

2 _ qi^ + (3lml _ 2 _ q^ + of^ 
1 — Pi 1 — «2 

Sp = {q+ + q-f = {qi + q2f ~ a2(3is - (qi + qa)^. 

Hereafter qi always stand for qi± etc. The conservation laws provides the relations 
among the introduced variables and thus the limits on ^2 variation 

pi = p- + p+, 02 = «- + «+, qi + q2 = q- + q+, (7) 

Sp + (qi + q2)^ n ^ , . 
< Pi ~ 0:2 < i. 

s 

As we shall see later the production amplitudes takes a simple form in terms of the 
following variables 

z±=P±/Pi, z- = l-z+, k± = q±-z±qi. (8) 

In what follows we will work in the kinematics of pionization region (similar approach 
was used in [|^), which give the main contribution to the total cross cross section of 
the process under consideration. For such kinematics the following limitations take 
place 

2 

— < aa ^ A < 1. 

s 

As a result the expression (^) can be represented as 

It is convenient to use the gauge invariance condition qi^^T^u = 0, which in view of (|^) 
entails 



^The spin of the colliding nuclei is unimportant because of the s-channel helicity conservation in 
the interaction of relativistic particles with Coulomb field. 



5 



UfJ, 



T — 



|qi 



qi 



A I3i ' ■ ' iqil 

Substituting this relation in (D and after a bit of Dirac algebra we find 



(10) 



^(1) _ |qil 



meiRi + iiQi + 22;+Qiei + 2|qi|z+z„_Ri 



P2 



(11) 



i?i = S(k_) - S{k+), Qi = k_5(k_) + k+5(k+), 5(k) 



k^ + 2+2_.qi^ + 



This representation for the amplitude is manifestly gauge invariant, i.e. its vanishes 
when qi — >• 0. Doing in ([TlD the replacement qi = q+ + q- — q2 it is easy to see that 
it vanishes also when q2 — >^ 0. 

For completeness we cite here the alternative expression for the amplitude which can 
be obtained if one use the gauge invariance constraint by another photon q2^T^i, = 0. 



R 



(1) 
(1) 



q2 r 

= mesi?! + eaQi + 2?/+Qie2 + 2\q2\y+y-Ri 

sa2 L 

R^ = S(l_) - S{\+ 



Pi, 



(12) 



S{1) 



1 



w? + q,2^y+y- 

The two factors 



0, Qi = 1+5(1h 

i± = q± - i/±q2, 



i_S(i_) 



y± = — , 

^2 



y+ + y- = 1- 



|qil 



|q2l 



qi^ + ml/Jf + mla^ 



are the familiar bremsstrahlung amplitudes for the photons with transverse momenta 
qi,2, which determines the luminosity of 77 collisions 0. The square of 'u(g_)i?|]^jf (g+) 
(up to the well known kinematics factors) gives the cross section of the process 7*7* — > 
e"'"e~. 

After a lengthy but straightforward calculation we obtain 



L(g_)/?|;]t;(g4 



= \sp{q- + m)i?|}J*(g+ - m)i?j|) 



(1) 



-z+z_q^ [{m^ + Az+z.q^)Rl + {zl + zi)Cl^ + Az+z.{z+ - 2_)i?iqiQi] . (13) 

The different terms in this expression have a transparent physical meaning in terms of 
the transverse and scalar (often called longitudinal) polarizations of the virtual photon 
7^ and there is one-to-one correspondence with the discussion of helicity structure func- 
tion in diffractive DIS and diffractive production of vector mesons [T^. Namely 
in the expression (|TBp the first term proportional to the corresponds to pair pro- 
duction by transverse photon with the sum of the lepton helicities A + A = ±1, the 
term (2;^ + 2;^)Qi^ provided the excitation of the pair with A -|- A = by the transverse 
photons. The term oc Azj^z_q^R\ in (|13D describes the pair production by scalar pho- 
tons, whereas the last term in (|13|) is the ST (LT) interference which vanishes upon 
the phase space integrations by 2+ and z_. 

Because of the exact conservation of the helicity of relativistic electrons in a Cou- 
lomb scattering, the above structure of the amplitude ([TT|) and its square ([T3|) will be 
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retained also in higher order. This important property will be widely used later on 
constructing the terms of higher orders in fine structure constant. 

The differential cross section of the process (|) is connected with the full amplitude 
by the relation 

d^ = ^El^l'dr. (14) 

The summation is carried over the final particles polarization. The phase space volume 
for four particles in the final state can be expressed through the variables introduced 
above and phase volume of the final leptons 

dr± = + (15) 

in the following way 

1 dVidV2dVdV.4. ^ , / ^ nfi\ 

= (2^ 2E[ 2^7 2^;^ + - - - - 

^^M^q,d^q.dr^. 



4(27r)8s(l-a2)(l-/3i) Pi 

In obtaining this expression we used the following relations for the final particles mo- 
menta ^ 

^ = 25{q^ - m^)d^q = sdadpd\5{sap - - m^). 
E 

The above expressions allows one to calculate the differential cross section of the process 
under consideration in the Born approximation with power accuracy in energy. 

The total cross section of pair production in the elementary process 7i72 — ^ e+e" 
drops rapidly as soon as one of the photons is off-mass shell, which imposes the upper 
limit of integration over the photons transfer momenta qi^,q2^ This limits are 

much smaller than the characteristic momentum where the formfactors are at work 

I ,2/ 2 6 10 0.25 ^ 

where Ra ~ 1.2^4^ fm is the nuclear radius, so that the effect of nuclear formfactors 
can be safely neglected. (It is not the case for the production of muon pairs or hadronic 
states though.) 

For this reason the evaluation of the total cross section in the case of electron pair 
production can be done in the point like nucleus approximation, i.e. using the cited 
above expressions. It can be shown that integrating the expression for the differential 
cross section ( p^ using the relations (^, (|13]) and (|16D one obtains the famous Racah 



formula for the total cross section 

a^ZlZl 28 



Tim? 27 



- 2.2L^ + 3.84L - 1.636 



(17) 



where L = ln(7i72). The total cross section of the process under consideration in 
the Born approximation is the growing function of the invariant energy (as the third 
power of logarithm in energy). This is the direct consequence of the fact that both 
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photons are quasi real, which bring to second power of logarithm in ([T^) and one power 
come from the common for all orders integration by rapidity variable jSi. Before going 
to calculation of properly Coulomb corrections let us do the following remark. The 
lepton pair created by two photons is produced in charge even i.e. spin-singlet state. 
The exchange of the odd number of photons between the Coulomb fields of heavy ions 
and produced leptons will lead to a pair creation in charge odd i.e. spin-triplet state, 
which doesn't interfere with Born amplitude. Nevertheless later on we have to take 
into account only the interference of the Born amplitude with the amplitude where the 
four photon exchanges take place between leptons and colliding particles (see Fig. ^). 
Now we can go to the main topic of present work, calculations of Coulomb corrections 
(CC). We begin from the simplest case which are provided by three photons exchange 
between the produced lepton pair and colliding particles. 



3 The Coulomb corrections. The spin-triplet 
(C-odd state) pair production 

We start with the discussion of the triplet pair production,which proceeds via the 
exchange of the odd number of photons between colliding particles and produced pair. 
The amplitude relevant to first CC of such type consist of two parts M(i2) = M^^^^+Mj-f^, 

where the amplitude ^(^2) corresponds to the case when the lepton pair is attached to 
colliding particle A by one photon and connects with B through two photons in all 

possible ways. The case when two photons are attached to A (one to B) is described by 

(2) 

the term M^^y As would be shown later this term can be obtained one from another 
by simple substitution, so it is enough to consider one of them. 

The amplitude Mj^-j is described by six Feynman diagrams which are depicted in 
the Fig. The only essential difference from the Born approximation is the additional 
photon exchange with four-momentum k = ap2 + (3pi + fc_L by which the integration 
has been done in the amplitude. To integrate over Sudakov variables a and (3 we use 
twelve FD instead of six which is relevant to this amplitude and introduce the statistical 
factor This trick permits us to provide the eikonal type integration over (3 using 
the identity 

J 2Txi \sl3 + iQ -s/3 + iO 

The convergence of the integration over a is provided by all six FDs obtained by 
permutation of all absorption points of exchanged photons by the pair. Closing the 
a integration contour to the upper half-plane one can see that only four of them are 
relevant. 
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Using the same approximation as in the Born case we obtain 



M, 



(1) 

(2) 



meiR2{k) + eiQ2{k) 



3? 



(1) 

(2) 



27r(qi^ 

|qil / d^k 1 



(19) 



S(3i J TT k2(q2-k)2 

+2z+Q2(k)ei + 2|qi|2;+2;_/22(fc)]P2 

|qil 



meiB!f} + eiQY'. + 2zj^Q,\^-.e-i + 2|qi|2;+z__R 



"-(2) 



'(2) 



(2)^ 



>(1) 
\2) 



P2, 



where in close analogy to amplitudes of diffractive DIS [|n 



i?2(k) 

Q2(k) = k_S(k_) 



5(k_) + 5(k+) - S(k_ - k) - 5(k+ - 

- k+S(k+) - (k_ - k)5(k_ - k) + (k^ 



k), 



(20) 



k)5(k+ -k). 



One can see that the spin structure of first CC is precisely the same as in the Born case. 



This is the direct result of lepton helicity conservation and as was shown in |17, 26 



is valid to all orders of perturbation series for the case, when one consider the single 
photon exchange between the pair and projectile A with any number of exchanges 
with target B and vice versa. Using the relations from Appendix A one can do the 
integration in (^) and obtains the analytical expression for 3^(2) • ■'■'^ general case this 
expression is rather cumbersome therefore we will present in the Appendix A only the 
result of integration in the case, when the single photon which is attached to projectile 
A is real (qi — > 0). The square of amplitude M^^^ can be obtained by analogy with 
the Born case 



2^+^-qi 



+ (4 + ^- 



(2) 



• (21) 



The amplitude M^^i^ ( two-photon exchange between the lepton pair and the nucleus 
Ai) is readily obtained from by following substitutions 



"2, 



^±^l/±, qi ^ q2, ei ^ e2, P2^Pi- 



(22) 



Let us estimate the considered CC contribution to the total cross section. This can be 
easily done if one restricted to the terms which grows with energy as the main power 
of the logarithm. In this case the CC behavior is defined by the second power of large 
logarithm L = ln(7i72) i.e. 



Pi 



d'qi 



qi 



IT 



Imin 



7i72 



so that to the leading logarithmical accuracy 

a'^ZfZliZf + Zi)L^ 



'odd 



(23) 



d^q_ 
TT 



d^q^ 



TT 



\2)) 



Q 



qi=0 



9 



Using the expressions for r\^^ 



(,^2) and QI2) from Appendix A we calculated this quantity 
with the result P ~ 3,6. We want to emphasize that the general expression (pjj) 
contains the leading logarithms contribution as well as the nonleading ones omitting 

2 

only terms which are suppressed by factor — and vanish in the high energy limit. 

The considered case qi ^ corresponds to WW approximation i.e. to the case 
when the scattered ion A is not detected and will scatter at very small angles. We see 
that for the case with semi-inclusive setup (when the scattered ion is fixed) we have 
much more complicated picture, which nevertheless can be described by cited above 
expressions. 

The contribution from the interference of the amplitudes M^^^) and M^^^j* is enhanced 
only by the first power of large logarithm L due to the boost effect and its impact 
on the total cross section will be consider further. Here we only cite this interference 
contribution to the charge odd part of the total cross section in the general form 



int 
odd 



a^ZfZlL 



d^qi d^q2 
vrqi^ irqi 



dsr 



(1) 

(2)- 



(24) 



4 The Coulomb corrections to the spin singlet 
(C-even) pair production 

The spin singlet pair production CC are provided by four photons exchange between the 
colliding ions and produced pair in all possible ways. Conventionally one can divided 
them in two set. The first one is determined by FD some of which are depicted on 



Fig. ^ a, b. The corresponding amplitude denotes by M13 = ^^''3) + M^^^^ represent the 
sum of all possible single attachments to one of the colliding ions with three photons 
exchanges between the lepton pair and another ion. The second set is determined by FD 

depicted on Fig. Rl c and include all double photon exchanges between lepton pair and 

(2) 

every colliding ions. This amplitude we write as M}J . Using the same approximations 



r{3) 



and technique as above for the amplitude M^^} we obtain 



(2)' 



M, 



(1) 

(3) 



IS 



{ATrayZiZl 



(3) 



1 



3! (27r)2(qi2 + /52m2- 

|qi 



(25) 



TT 



TT ki2k2^(q2 — ki — k2 



X 



meiR-iik) + eiQ^lk) + 2z+Q3(A;)ei + 2\qi\z+z_Rs{k) 



P2 



|qi| 



meii^lgj + eiQ(3) + 2z+Q|3|ei| + 2\q-i_\z+z_R^^ 



P2, 



where in close analogy to amplitudes of diffractive DIS 

i?3(ki, k2) = -S(k_) + S(k_ - ki) + S{k_ - k2) + S(k_ - q2 + ki + k2) 



(26) 



Q3(ki,k2 



k_5(k_) - 

- (k_ - q2 



ki 



ki)5(k_ - ki) - (k_ - k2)5(k_ - k2) 

h k2)^(k_ - q2 + ki + k2) - [k_ ^ k+ 



One can see that as in the case of third order CC considered above the fourth order 
amplitude of first kind (i.e. with one single attachment) has the same spin structure as 
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the Born one and can be obtained from it by simple replacement (-Ri, Qi) ^ (i?3, Q3). 
To get the amplitude M^-^^ it is enough to do the replacement (^) in the above expres- 
sions. It is easy to check that the amplitude M12 vanish when any of photons momenta 
qi or q2 goes to zero, so that the gauge invariance is provided. Finally, in the WW limit 
qi,q2 — ^ one obtains from (125|), (|26D the expression which coincide with respective 



one from p7 |. 

The main contribution to the total cross section come from the interference of 
the M13 with the Born amplitude M^^y To the leading logarithmic accuracy this 
interference can be read 



^(1) + ^(3) 



Pi, 



(27) 



Pi 



d^q2 d^q_ 



TT 



qi=0 



Compare this expression with the cross section in the case of C-odd CC ( p3| ) one can 
see that the both terms are of the same order in fine structure constant and have the 
same energy dependence. Moreover we make a point that Pi = P. Indeed, the product 
RiR^^^ contains 16 terms and making the proper shifts of the integration variables. 



'(3) 

one can establish one-to-one correspondence to the 16 terms in ( R 



^2) 



Similarly the 



same can be done with the vector products QiQ^gj and ( Q(2j ) • This can be regarded 



as a manifestation of the well known AGK rules ||22[ . Really, the two photon exchange 
production of lepton pair can be viewed as a diffraction excitation of the photon, 
whereas the 3-photon exchange can be viewed as an absorption correction to the one- 
photon exchange. Furthermore, the total cross section for pair photoproduction in the 
Coulomb field reads fSBI 



28Z^a ^ 



109 
I2" 



^ nin^ + iZaf) 

n=l ^ \ / / 



i2t 



which for the relevant order gives 

P = Pi = ^C(3) ^ 3.6, C(3) = f2-3 = 



n=l 



As was mentioned above this result was verified by numerical calculations. 

Returning to the second set of the fourth order CC (two photons exchange between 
the pair and every heavy ion) we calculated all respective terms in accordance with the 



FD (see Fig. y c) and obtained the matrix element M}J in the following form 



(2) 



3? 



(2) 

(2) 



M, 

d^kl 



(2) 
(2) 



TT 



isa^2\ZiZ27iyu{q_)^[i\v{q+), 

I ^ n , p \d(2) 



(29) 



d^k2 



The structure 3fJ^2) so called up-down symmetry, which is decoded in the permu- 
tation operator Pud, which acts as follows 

PudP{pi,P2,a,P,qi,ki,q2,k2) = F{p2,pi, P,a,q2,k2,qi,ki). (30) 
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The quantity Rj^^!^ is the sum of all terms corresponding to all possible FD for this case 

-R(2) = -R1243 + -R1324 + -R243I + -R3124 + -R1342 + -Rl423- (31) 

The structures Rigki are cited in Appendix B. Unfortunately we do not succeed in 
trying to write down the result for spin structures in this case as it has be done above 
for the terms where there is a single attachment of exchanged photon to one of the 
colliding ions. 

Nevertheless we cite here the leading contribution to the C-even state pair creation 
cross section which come from the interference of this amplitude and the Born term 

<L = ^^^^ J / "-''I + - '''^ 

This expression is the same order in fine structure constant as the considered above 
interference between the Born amplitude and first set of fourth order term (see (p7D), 
but it has a one power less in its energy dependence. 

On the other side the first power dependence on L and the same order in a has 
contribution to the total cross section arising from the interference between the ampli- 
tudes -^(^2) ^^"^ -^(1) (^^^ (iH)) which corresponds to pair production in C-odd state. 
From very general arguments this contributions have to cancel each other within the 
logarithmical accuracy. This statement^ follows from the fact that the relevant (sixth 
order in a) elastic amplitude for AB scattering at zero angle (see Fig. ^ has only single 
logarithmical enhancement factor due to the boost freedom of a light-light scattering 
block, whereas the s-channel discontinuity of this amplitude has no logarithmical en- 
hancement at all. Therefore in the total cross section of pair production the terms 
proportional to (ZiZ2a^)^L are absent. It is interesting to check this statement by 
straightforward calculations, which will be done elsewhere. Moreover it seems very 
natural that the contributions of the form {Z\a)^^^^{Z2Cif'^^^ are absent in the total 
cross section in any orders of perturbation theory. This statement is in the complete 
agreement with the results recently obtained for the Coulomb corrections ||28|| . 

5 Concluding remarks 

In this work we have considered lepton pair production in the Coulomb field of two 
highly relativistic heavy ions. Using the powerful technique of Sudakov variables,which 
allows one to obtain the results which are correct with the power accuracy in energy, 
we get the first terms in perturbation expansion of amplitude (up to the fourth order in 
fine structure constant) and show that they are finite and gauge invariant. We obtain 
the simple analytical expressions for the Born amplitude and relevant cross section, 
which allows one to calculate the lepton pair yield in the wide interval of kinematic 
variables. It is shown that the terms in the amplitude, which corresponds to FD with 
at least one single exchange between created pair and one of the ion, have the same 
spin structure in all orders in a and we proposed the simple rules to build them. 

In every order of perturbation theory we analyze the obtained CC isolating the 
leading in energy terms and show that remarkable cancellations among the different 
terms in the total cross section take place. Finally the numerical estimates of the main 
contribution from the CC in the total cross section have been done. 

^We are grateful to L. Lipatov for discussion on this issue. 
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Appendix A 

We give here the result of two-dimensional integration by loop momenta for scalar and 
vector structures determining the first CC 



/(q_ 



d^k 



TT (k2 + A2)((q-k)2 + A2)((q_ -k)2 + m2 



1 



q 

+ Ai In 



q 



+ 



qn 



+ ^2 In 



A2 A L(g-) 



- 1 



In- 



m 



(A.l) 



^1 

^2 



qq- 
qq+ 



1 

-r + 



1 



+ 



+ 



A 







^^qq+ 


{«-> 


(9+) 




H 




m2qq_ 


{«-> 


(9+) 





(A.2) 



i(q 



d^k 



TT (k2 + A2)((q-k)2 + A2)((q_ -k)2 + m2) 

1 



q 1 (?+) 



q2(g+) 
1 

-q 



A 



q- 



q- 



(A.3) 



where we use the notation (a) = a? + m? and besides 

A = -m2q2, q = q+ + q_. 

Using these integrals and the expression (^) we obtain for scalar combination which 



enters M, 



(1) 

(2) 



R 



(1) 

(2) 



d^k R2{k) 



TT 



2 
A 



k2(q - k)2 



q 



(A.4) 



R2{k) 



1 



^ 

{<!-) (?+) 
1 



In 



m2q2 



q^ v(?+) 



In 



1 



{k-q^) {k - q+)' 

For the case of the vector structures we obtain 
d^k Q2(/c) 



^(2) 



TT k2(q-k)2 

,2 / 1 



q- -q 
q 



Q2(A;) 



+ 

q 

+ 2 

q- 



A \{q-) 



In 



m2q2 q2 \{q^) 



In 



(A.5) 



m2q2 



In- 



q- 



In- 



-) 



+ 



m2 m2 
k — q q+ k — q 



{k-q-) {k-q+y 
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Appendix B 



We cite here the spin structures which determines the fourth order amphtude of second 
class (see expression (|3^)) 

1 Pijq- - qi + m)p2 

-^1243 = --^ — ■ ■ ■ [D-i] 



-Rl324 

-R243I 
-R3124 
-R1423 



R 



2413 



S |±(g_) + (g_-g,) 

1 Pi{q- - ki + m)p2{q- - h - k2 + m)pi{-q+ + q2 - k2 + m)p2 
^{q-){q--ki-k2) + {q--ki){-q+ + q2-k2) 

1 pi{q- - qi + ki + m)p2{-q+ + ki + m)pi 
a-{-q+ + ki) + a+{q^ - qi + ki) 

1 P2{q- ~ h + m)pi{-q+q2 - h + m)p2 
/3+{-q^+k2)+/3-{q^-qi-k2) ' 
1 Pi{q~ - ki + ■m)p2{-q+ + qi + k2 - ki + m)pi{-q+ + k2 + m)p2 
^{q-){-q+ + qi + ki-k2) + {q--k^){-q+ + k2) 

1 Piiq- - qi + ki + 'm)p2{-q+ + ki + k2 + m)pi{-q_^ + h + m)p2 
^{q^){-q+ + k2 + k^) + {-q+ + k2){q2-q+ + k^) 



The symmetry property is caused by the charge-conjugation symmetry of the amph- 

/2) 

tude. There are 24 FD contributing to M^^y Instead of them it is convenient to 
consider 24 * 2 * 2 = 96 FD which take as well the permutations of emission and ab- 
sorption points of exchanged photons to the nuclei. The result must be divided by 
(2!)^. This trick provides the convergence of ai,/32 integrals. All the 24 FD, describing 
the interaction of 4 virtual photons with the pair components are relevant providing 
the convergence of 02 integrations. When closing the contour of integration say in 
upper planes only 12 from them become relevant. 



16 




Pi p'l 




P2 P2 P2 P2 

(a) (b) 

Fig. 1: The lowest order Feynman diagrams for the two-photon pair production in the 
process A + B—>-A + B + e+e". 





(d) 



(e) 



(f) 



Fig. 2: Gauge invariant set of Feynman diagrams for spin-triplet production which 



determines the matrix element M, 



(2) 



17 



Pi 



p'l 






P2 



P'2 



(a) 



(b) 



(c) 



Fig. 3: Some Feynman diagrams for the four-photon pair production in the process 
A + B ^ A + B + e+e-. 



Pi 



A 



B 



P2 



(a) 



p'l 





P2 



(b) 





(c) 



(d) 



Fig. 4: Different s-chanel cuts for forward elastic scattering in sixth order of perturba- 
tion theory. 
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